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ABSTRACT 

The size distribution in the cold classical Kuiper belt can be approximated 
by two idealized power laws: one with steep slope for radii R > R* and one 
with shallow slope for R < R*, where R* ~ 25-50 km. Previous works suggested 
that the SFD roll-over at R* can be the result of extensive collisional grinding 
in the Kuiper belt that led to the catastrophic disruption of most bodies with 
R < R*. Here we use a new code to test the effect of collisions in the Kuiper 
belt. We find that the observed roll-over could indeed be explained by collisional 
grinding provided that the initial mass in large bodies was much larger than the 
one in the present Kuiper belt, and was dynamically depleted. In addition to the 
size distribution changes, our code also tracks the effects of collisions on binary 
systems. We find that it is generally easier to dissolve wide binary systems, 
such as the ones existing in the cold Kuiper belt today, than to catastrophically 
disrupt objects with i? ~ i?*. Thus, the binary survival sets important limits on 
the extent of collisional grinding in the Kuiper belt. We find that the extensive 
collisional grinding required to produce the SFD roll-over at R* would imply 
a strong gradient of the binary fraction with R and separation, because it is 
generally easier to dissolve binaries with small components and/or those with 
wide orbits. The expected binary fraction for R < R* is <0.1. The present 
observational data do not show such a gradient. Instead, they suggest a large 
binary fraction of ~0.4 for R = 30-40 km. This may indicate that the roll-over 
was not produced by disruptive collisions, but is instead a fossil remnant of the 
KBO formation process. 

Subject headings: Kuiper belt: general 
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1. Introduction 

The cold Classical Kuiper Belt, hereafter cold CKB, is a population of trans- Neptunian 
bodies dynamically defined as having orbits with semimajor axis a — 42-48 AU, perihelion 
distances that are large enough to avoid close encounters to Neptune, and low inclinations 

Cold CKB objects (cold CKBOs) show several properties that distinguish them from 
other populations in the trans-Neptunian region. Specifically, cold CKBOs have distinctly 
red colors (Tcgler & Romanishin 2000) that may have resulted from space weathering of 
surface ices, such as methanol (e.g., Schaller 2010), that are stable beyond 30 AU. A large 
fraction of 100-km-class cold CKBOs are binaries with wide separations and similar size 
components (Noll ct al. 2008a,b). These binaries are practically absent in the dynamically 
hot, resonant and scattered populations. The albedos of cold CKBOs are generally higher 
than those of the dynamically hot CKBOs (Brucker et al. 2009). And finally, the magni- 
tude distribution of cold CKBOs is markedly different from those of the hot and scattered 
populations, in that it shows a steep slope for the large objects: E(mR)dmR ~ io°"*R with 
a ~ 0.82, where S(mR) is the number of objects per square degree with R magnitude 
(Eraser et al. 2010). 

The magnitude distribution of CKBOs shows a roll-over to a shallower slope at mag- 
nitudes ~ 25 (Bernstein et al. 2004, sec Petit et al. 2008 for a review). This feature 
is thought to occur due to a genuine change in the slope of the size frequency distribution 
(SFD), because extreme albedo variations would need to be invoked if the SFD slope were 
constant over m^. Using albedo pv = 0.2 for the cold CKBOs (Brucker et al. 2009) and 

— 25 (Fuentes et al. 2009), the roll-over radius is R* = 37 km. Allowing for uncertainties 
in Pv and m^, it is probably safe to assume that 25 < i?* < 50 km. 

For rriR < m^, a ~ 0.82 from Fraser et al. (2010) implies N{R)dR ~ R'^dR with 
g = g> ~ 5 for > i?* (q — 5a + l). For tbr > m^, a ~ 0.2-0.3 from Fuentes & Holman 
(2008) and Fuentes et al. (2009) imphes that q = 2-3 for R < R*. These SFD 

slope estimates will be used to constrain our model (§3). Note that the shallow SFD slope 
for mR > is an empirical approximation that fits observations only up to the limiting 
magnitude tbr ~ 27 (Fuentes et al. 2009) or equivalently down to i? ~ 15-20 km, depending 
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on albedo. The actual SFD can be wavy, just like the SFD of the asteroid belt at these radii 
(e.g., Bottke et al. 2005). Indeed, the occultation data from Schhchting et al. (2009) may 
indicate that the number of i? > 250 m KBOs is significantly larger than what would be 
expected from extrapolating the SFD with g < 3 from R* down to R = 250 m. 

The SFD roll-over at R* may be telling us something important about the formation and 
evolution of cold CKB. As first suggested by Pan & Sari (2005, hereafter PS05), the shallow 
SFD slope for R < R* can be the result of extensive coUisional grinding in the cold CKB. 
Indeed, the SFD slope index below i?*, g< ~ 2-3, is similar to the slope index expected for a 
collisionally evolved population that reaches the equilibrium slope (Dohnanyi 1969, O'Brien 
& Greenberg 2003) t 

PS05 carried out order-of-magnitude analytic calculations to support their argument. 
They used an idealized two-slope SFD with g< ^ 3 and g> = 5, and estimated R* by 
postulating that bodies of radii R < R* are disrupted over 4.5 Gy. As presented, their 
results were meant to apply to the case in which the number of bodies with R > R* has 
not changed over 4.5 Gy, and was equal to the number of bodies in the present KB (but see 
below). They found that R* = 20-50 km, in agreement with observations, with the exact 
value of R* mainly depending on the assumed disruption law. 

We were unable to reproduce these analytical results with the numerical code described 
in §2. Specifically, using PS05's assumptions, we always obtained the break well below 
10 km (see §3.1). Since our code was thoroughly tested, and gives results that are essentially 
identical to those of other numerical codes (Weidenschilling et al. 1997, Stern & Colwell 
1997, Kenyon & Bromley 2001, Bottke et al. 2005, Fraser 2009), one might wonder why it 
fails so badly in reproducing the PS05 results. While some unspecified numerical problems, 
common to all tested codes, can be responsible, there are also a few issues with the PS05 
calculations that need to be clarified. 

PS05 claim that they normalize the population at R* to be equivalent to KB with 4 x 10^ 
bodies larger than 100 km. A thorough reading of the description of their normalization 



^Dohnanyi's equilibrium slope q — 3.5 (Dohnanyi 1969) applies to a situation in which the object's 
strength is independent of its size. In the gravity regime of impacts, applicable the size range relevant here, 
the strength increases with size and q ^ 3 (O'Brien & Greenberg 2003). 
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procedure reveals, however, that they mean R > 100 km rather than diameter D > 100 km. 
Using this normahzation, they therefore have 4 x 10^ bodies with R > 100 km, which imphes 
~6.4 X 10^ bodies with i? > 50 km, if R* < 50 km and g> = 5. For comparison, observations 
indicate that there should only be 30,000 to 50,000 KBOs with i? > 50 km (Jewitt & Luu 
1995, Trujillo et al. 2001, Fuentes & Holman 2008). PS05 therefore effectively assume a 
population of R > R* bodies that is a factor of ~16 larger than the one existing in the 
KB today. Since, according to PS05's eq. (8), R* scales with f^^^, where / = 16 is the 
augmentation factor, we find that R* should be about a factor of 2.5 lower than what was 
claimed in PS05, or R* = 8-20 km, if / = 1. 

The second issue is related to the PS05's disruption laws. To obtain large values of 
R*, they assumed extremely weak disruption laws {(3 = 1.5 in the notation of their eq. 5), 
where the specific impact energy needed to catastrophically disrupt and disperse an object 
of size R, Q}j{R), was up to ~ 10'^ times lower than the one estimated for water ice by Benz 
& Asphaug (1999) (Fig. [T]). Such extremely weak disruption laws are probably unrealistic. 
More realistic disruption laws, such as those advocated by Stewart & Leinhardt (2009) and 
Leinhardt & Stewart (2009), are similar to the strongest disruption law considered by PS05 
(/3 = 3 in their eq. 5). With these laws, R* < 10 km according to PS05 (for / = 1), which 
is more in line with the results of our numerical code. We therefore find that either the 
number of large objects was much larger than the one in the present KB, implying faster 
initial grinding, or the observed roll-over at R* = 25-50 km must be due to something else. 

The SFD roll-over can be a signature of the accretional processes that were active 
during the KBO formation (see Kenyon et al. 2008 and Chiang & Youdin 2010 for reviews). 
Interestingly, the asteroid belt, thought to have formed by similar processes, also shows a 
roll-over at i? ~ 50 km. Given better constraints that we have on the fragmentation processes 
in the asteroid belt (e.g., survival of the Vesta's crust, number of asteroid families, etc.), it 
has been established that the roll-over was not produced by coUisional grinding (Bottke et 
al. 2005). Instead, the asteroid SFD roll-over was likely created by accretional processes, 
and constrains them in important ways (Morbidelli et al. 2009). The distinction between 
fragmentation and accretion signatures is more difficult to make in the Kuiper belt, where 
we are lacking good constraints. 

The binary KBOs can provide an interesting constraint on the amount of collisional 
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grinding in the Kuiper belt. Recent observations indicate that ~30% of 100-km-class cold 
CKBOs are binaries (Noll et al. 2008a,b; >0.06 arcsec separation, <2 mag magnitude 
contrast). The properties of known binary KBOs differ markedly from those of the main- 
belt and near-Earth asteroid binaries (Merline et al. 2002, Noll et al. 2008a). The 100-km- 
class binary KBOs identified so far are widely separated and their components are similar 
in size. These properties defy standard ideas about processes of binary formation involving 
coUisional and rotational disruption, debris re-accretion, and tidal evolution of sateUite orbits 
(Stevenson et al. 1986). They suggest that most binary KBOs are remnants from the earliest 
days of the Solar System. Indeed, all models developed so far postulate that binary formation 
was contemporary to the formation of KBOs themselves (Weidenschilling 2002, Goldreich et 
al. 2002, Funato et al. 2004, Astakhov et al. 2005, Nesvorny et al. 2010). 

The survival of binary KBOs after their formation is an open problem. Petit & Mousis 
(2004) estimated that several known binary KBOs, such as 1998 WW31, 2001 QW322 and 
2000 CF105, have hfetimes against coUisional unbinding that are much shorter than the age of 
the solar system. These estimates were based on the steep SFD adopted by Petit & Mousis 
that was extended down to R — 5 km. This assumption favors binary disruption, because 
of the large number of available impactors. When we update Petit & Mousis' work with the 
more recent estimates of the SFD in the Kuiper belt (e.g., R* = 25-50 km for cold CKBOs; 
see discussion above), we find that a typical 100-km-class wide binary CKBO is unlikely to 
be disrupted over the age of the solar system (^1% probability), except if the KB (or its 
source population) was more massive/erosive in the past. 

Here we use the observed binary fraction in the cold CKB to determine how it limits 
the amount of coUisional grinding in the Kuiper belt. We make different assumptions on the 
initial state and history of the cold CKB, and the related populations, and identify cases 
that lead to the SFD break of cold CKBOs at R* — 25-50 km. We then evaluate the survival 
of binary CKBOs in each of these scenarios. If the roll-over at R* — 25-50 km was produced 
by coUisional grinding, we find that the binary fraction should show a strong gradient with 
radius and binary separation, because it is generally easier to dissolve binaries with small 
components or those with very wide orbits (§3). The absence of such a gradient would 
indicate that the roll-over is accretional. 
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2. Modeling CoUisional Evolution and Binary Survival 
2.1. CoUisional Evolution Code 

Our collisional modeling simulations employ Boulder, a new code capable of simulat- 
ing the collisional fragmentation of multiple planetesimal populations using a statistical 
particle-in-the-box approach. It was constructed along the lines of other published codes 
(Weidenschilling et al. 1997, Kenyon & Bromley 2001). A full description of the Boulder 
code, how it was tested, and its application to both accretion and collisional evolution of 
the early asteroid belt, can be found in Morbidelli et al. (2009). Examples of its previous 
use for the asteroid belt, Hildas, Trojans, irregular satellites, and primordial trans-planetary 
disk are described in Levison et al. (2009) and Bottke et al. (2010). 

The code's procedure for modeling impacts is as follows. For a given impact between a 
projectile and a target object, the code computes the specific impact energy Q, defined as 
the kinetic energy of the projectile divided by the target mass, and the critical impact energy 
Q*j^, defined as the energy per unit target mass needed to disrupt and disperse 50% of the 
target (e.g., Davis et al. 2002). For reference, Q < values correspond to cratering events, 
Q ^ Q*j^ correspond to barely catastrophic disruption events, and Q > Q}) correspond to 
super-catastrophic disruption events. 

For each coUision identified by the code, the mass of the largest remnant is computed 
from the scaling laws found in hydrodynamic simulations of impacts (Benz & Asphaug 1999, 
Leinhardt & Stewart 2009, Stewart & Leinhardt 2009). The mass of the largest fragment 
and the slope of the power-law SFD produced by each coUision is set as a function of Q/Q^ 
by empirical fits to the hydrocode results of Durda et al. (2004, 2007) and Nesvorny et 
al. (2006). See Bottke et al. (2010) for the explicit definition of these fits. These results 
apply to monolithic target bodies. We also tested approximate scaling laws for impacts on 
pre-fragmented and rubble-pile targets. Again, these scaling relations were drawn from the 
fits to hydrocode impact simulations (Benavidez et al. 2011). 

The Q}) function was assumed to split the difference between the impact experiments 
of Benz & Asphaug (1999), who used a strong formulation for ice, and those of Leinhardt 
& Stewart (2009), who used the finite volume shock physics code to perform simulations 
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into what they describe as weak ice. To do this, we divide the Benz & Asphaug (1999) 
strong-ice function by a factor, fq. We typically used fg = 3, 5, and 10. Note that 
because we sampled a broad section of parameter space, we chose not to include still more 
complicating factors (e.g., QJ-, may vary with impact velocity, etc.). The bulk density was 
set to p = 1 g cm^^. 

The main input parameters for the Boulder code are: the (i) initial SFD of the simulated 
populations (see §3); (ii) intrinsic collision probability. Pi, defined as the probability that 
a single member of the impacting population will hit a unit area of a body in the target 
population over a unit of time, and (iii) mean impact speed, Vi. For collisions between 
present-day CKBOs, we used Pj = 4 x 10~^^ km~^ yr~^ and Vi = 1 km s~^ (Davis & 
Farinella 1997, Dell'Oro et al. 2001) @ See §3 for our specific choices of Pj and Vi for other 
populations. 



2.2. Binary Survival Code 

Petit & Mousis (2004) identified three main processes that can dissolve KBO binaries: 
(1) One of the component is hit by a small impactor. While the component survives essen- 
tially intact, except for a new crater on its surface, the linear momentum transferred from 
the impactor imparts a 'kick' on the velocity vector of the binary orbit. If the kick is large 
enough, the component becomes unbound from its companion. (2) One of the components 
can be shattered by a large impactor. The fragments produced the object's breakup are ex- 
pected to escape on unbound trajectories because the ejection speeds (~100 m s~^) largely 
exceed that of the binary orbit (~1 m s~^). (3) The binary system has a close gravita- 
tional encounter with another KBO. The tidal gravity of that object can unbind the binary 
provided that the object is massive enough (e.g.. Stern et al. 2003). 

Petit & Mousis (2004) found that mechanism (1) is by far the most efficient way to 
dissolve binaries in the present KB. Thus, we focus on modeling (1) in this work. Mechanism 



^Note that, according to PS05, Pi can be approximated by ~ ^/A, where ft is the typical orbital angular 
velocity in the Kuiper belt, and A is the Kuiper belt's area in the plane of the solar system. With il = 
0.022 yr~^ and A = 1200 AU^, PS05 thus effectively have Pi = 10"^^ km"^ yr"\ a value larger by a factor 
of 2.5 than the one adopted in this work. 
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(2) is also considered, but we find that its effects are negligible compared to (1), except if 
extremely weak (and probably unrealistic) disruption laws are adopted. This result stems 
from the fact that it is generally easier to dissolve a wide KBO binary orbit by an impact than 
to physically disrupt a 100-km object. Mechanism (3) could have been important during 
the early phases of the Kuiper belt evolution provided that the cold CKB overlapped with 
a population of numerous, very massive KBOs. We do not model (3) in this work because 
we do not yet have a detailed understanding of these early stages. 

We will assume in the following that small impacts with Q <^ Q}^ can be treated as 
inelastic collisions. Thus, we will ignore any linear momentum that can be potentially carried 
away by ejected fragments. In this approximation, and assuming that <^ m, where 
and m are the impactor and binary component masses, the velocity vector of the binary 
orbit, V, will change by 



See DeirOro & Cellino (2007) for a discussion of the linear momentum transfer for different 
impact angles, and in the case where the escaping ejecta affect the linear momentum budget. 

The binary components are assumed to have equal mass, which should be a good ap- 
proximation for the real binaries in the cold CKB (see §1). This assumption is conservative 
in the sense that it is generally harder to disrupt a binary with equal-size components than 
the one in which the secondary is smaller than the primary (if other parameters are the 
same). By using equal-size binaries we may thus slightly underestimate the real decay rate. 

The radial, tangential and normal components of 5v will be denoted by Svr, 6vt and 
Svz, respectively, in the following. With this notation, the semimajor axis a, eccentricity e 
and inclination / of the binary orbit will change by 




(1) 



Sa 



— [Svt (1 -I- e cos /) + 6vji e sin /] , 



(2) 



<R sin / + Svt (cos f+ ^'^ ) , 

\ 1 + e cos / / _ 



(3) 



S{cosI) 



rj sin(w + /) 

— ovz — 7 , 

na I + e cos / 



(4) 
(5) 



where u is argument of pericenter, / is the true anomaly at the time of impact, and rj — 
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Vl - e2 (e.g., Bertotti et al. 2003). 

Our statistical code does not deal with the detailed geometry of each individual im- 
pact. Instead, it follows the mean quadratic changes of orbital elements produced by the 
average effect of impacts with different orientations relative to the binary orbit, and times 
corresponding to different orbital phases. Specifically, we assume that the impact direction 
is isotropic and impacts are not correlated with the binary orbit phase, which should be the 
case for the well-mixed system that we deal with. The algorithm is as follows. 

Assuming that Vj = Vi n, where Vi is the characteristic impact speed and n is the unit 
vector with isotropic orientation, we have 

{{Svnr)-{{5vrr}-{{Sv,r)-l^v^, (6) 

and 

{5vr 5vt) = {Svr 5vz) = {Svt Svz) = , (7) 

where the average was taken over n's orientations. Speed Vi can be determined as the r.m.s. 
value of impact speeds with the distribution that is characteristic to the studied population 
(e.g., Davis & Farinella 1997). 

In the next step, we average over the phase of the binary orbit at which the impact 
occurs. We obtain 

((5a) ) 4 frriiV. 



2 

"t '-'I 



0? 3 V mv 
5 



(8) 



. . (9) 
^ ' 6 \ mv / 

and 



\2rf \ mv 

Since the orbital speed of a KBO binary, v = na, where n is mean motion, is typically 
~1 m s~^, while ~ 1 km s~^, the binary orbit can be dissolved by a single impact 
generating Sa/a, Se or S{cosI) of the order of unity, if mi/m > v/vi ~ 10~^. 

In addition, the cumulative effect of impactors with mi < 10~^m can also be important 
as it leads to a random walk in orbital elements. It can be shown that 



{5a Se) = {Sa S (cos /)) = {Se S (cos /)) = 



(11) 
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Thus, the changes in different orbital elements due to small impacts are uncorrelated and 
can be treated separately. 

The binary system is assumed to become dissolved either if e > 1 or if the semimajor 
axis exceeds some critical value, Ocrit- We set Ocrit = 0.5i?Hiii as a rough limit dictated by the 
Hill stability criterion (e.g., Donnison 2010) and numerical integrations (e.g., Nesvorny et al. 
2003). By experimenting with the code we find that the main channel of binary disruption 
is reaching a > acrit in one or a small number of collisions. Accordingly, the results are 
insensitive to the initial distribution of e. 

Note that the inclination changes cannot result in binary splitting directly, but, if cou- 
pled to effects arising from the Kozai dynamics (Kozai 1962), they can also be important. 
We do not consider the coupling of coUisional effects with dynamics of inclined binary orbits 
in this work. Therefore, the initial direction of the binary angular momenta does not need 
to be specified. 

2.3. Boulder with Binary Module 

The binary module was inserted in the Boulder code. This was done by attaching 
additional data structures to each size bin. These data structures describe the initial distri- 
butions of a, e and i of binary orbits, and track how these distributions change over time. 
As the population of impactors evolves with time due to collisional fragmentation, the code 
calculates the time-dependent rate of change of binary orbits, and evolves them according 
toEqs. (IH1)-([10]). 

3. Results 
3.1. PS05 Case 

To start with, we illustrate the case studied by PS05. In this case, the KB is assumed 
to evolve in isolation over r = 4.5 Gy. The initial population is given a two-slope SFD 
with g< = 3, g> = 5, and R* between 1 km and 50 km. The total number of objects 
with i? > 50 km is normalized to A^o = 50, 000, which is roughly the estimated number of 
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objects in the present KB. We also consider cases with = JmNq, where /tv is the scahng 
parameter. For example, assuming pv = 0.2, p = 1 g cm~^, g< = 3, g> = 5 and R* = 37 
km, ~ 0.4 gives the total mass ~ O.OlMg, where Me is the Earth mass, which is what 
the observations seem to indicate for the cold CKB (Fuentes & Holman 2008), although this 
value is still quite uncertain. On the other hand, Jn = 13 corresponds to the population 
used in PS05 (see discussion in §1). 

Figure H] shows the result of collisional grinding for fi^ = l and two different values 
of /q- /q = 3, roughly corresponding to the PCOS's strong disruption law, and /q = 500, 
roughly corresponding to the PCOS's weak disruption law (see Fig. [1]). We used Pi = 
4 X 10~^^ km~^ yr^^ and Vi = 1 km s~^. The initial distribution was set so that g< = 3 for 
R < 1 km and g> = 5 for i? > 1 km. For p = 1 g cm~^ this gives the total initial mass of 
0.84 Me. After having evolved the population with the Boulder code over r = 4.5 Gy, the 
remaining masses were 0.056 Me for /q = 3 (Fig. [2|i) and 0.027 Me for fq = 500 (Fig. [2b). 

For /q = 3, the SFD slope just below ~50 km becomes slightly shallower over time, 
mainly due to the effects of large cratering impacts, and reaches g ~ 4 at r = 4.5 Gy. This 
final slope is significantly steeper than the present slope of cold CKBOs at these radii. With 
fq = 500, a sharp roll-over to a very shallow slope develops at i? ~ 10 km, because most 
objects with i? < 10 km suffer catastrophic disruptions. Thus, even with the unrealistically 
weak disruption law, the SFD break is still significantly below the actual roll-over radius in 
the present cold CKB. This result is insensitive to the specific choice of model parameters 
related to the generation of fragments, resolution, and plausible changes of Pi and/or fj. 

We therefore find that the observed break at R* = 25-50 km cannot be produced by 
collisional grinding, except if the number of objects with R > R* was much larger in the 
past, and was depleted by dynamical processes, or if the cold CKBO overlapped with a much 
larger population of impactors in the past. We consider these possibilities in the following 
text. 



3.2. Adopted Model 



Levison et al. (2008a) proposed that most of the complex orbital structure seen in the 
KB region today (see, e.g., Gladman et al. 2008) can be explained if bodies native to 15- 
35 AU were scattered to >35 AU by eccentric Neptune (Tsiganis et al. 2005). If these outer 
solar system events coincided in time with the Late Heavy Bombardment (LHB) in the inner 
solar system (Gomes et al. 2005), binaries populating the original planetesimal disk at 15- 
35 AU would have to withstand ~600 My before being scattered into the Kuiper belt. Even 
though their survival during this epoch is difficult to evaluate, due to major uncertainties in 
the disk's mass, SFD and radial profile, the near-absence of wide and equal-sized binaries 
among 100-km-sized hot classical KBOs (Noll et al. 2008a,b) seems to indicate that the 
unbinding collisions and scattering events must have been rather damaging. 

The contrasting characteristics of cold CKBOs discussed in §1 may indicate that the cold 
CKBOs formed in a relatively quiescent environment at >40 AU rather than having been 
scattered to their current orbits from <35 AU, because more resemblance between different 
trans- Neptunian populations would be expected in the latter model. The in-situ formation 
of cold CKBOs is also supported by the results of Parker & Kavelaars (2010) who showed 
that some of the widest binaries observed in the cold CKB today would probably not survive 
scattering encounters with Neptune. 

To understand the collisional history of cold CKB, we consider both the prc-LHB and 
post-LHB epochs. For the pre-LHB epoch, we assume that the cold CKBOs evolved in rela- 
tive isolation at 40-50 AU, where most collisions occurred between cold CKBOs themselves. 
The results of these pre-LHB simulations should also apply, with minor modifications, if the 
cold CKBOs formed closer in, as far as the source population of cold CKBOs can be consid- 
ered as a closed collisional system. The coupling of CKBOs to the scattered trans-planetary 
disk during the LHB is not considered here because Benavidez & Campo Bagatin (2009) 
showed that collisional fragmentation during this stage was unlikely to produce substantial 
changes in the SFD. 
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3.3. Post-LHB Epoch 

We start by discussing the collisional evolution of the Kuiper belt after LHB. We consider 
collisions between two CKBOs, one CKBO and one scattered disk object, and two scattered 
disk objects. The scattered disk is massive initially and becomes dynamically depleted over 
time, with an estimated depletion factor of 100-250 over 4 Gy (Levison & Duncan 1997, 
Dones et al. 2004, Tsiganis et al. 2005). This dynamical depletion is taken into account in 
Boulder by gradually decreasing the number of scattered disk objects in all size bins. 

The collisional probabilities, impact speeds, initial SFDs and dynamical decay rates 
were taken from Levison et al. (2008b). Specifically, we used Pj = 4 x 10~^^ km~^ yr^^ and 
Vi = 1 km s~^ for CKB collisions, Pj = 2 x 10"^^ km"^ yr~^ and Vi = 1.5 km s~^ for CKB- 
scattered-disk collisions, and Pj = 10~^^ km^^ yr^^ and f j = 3 km s~^ for collisions between 
scattered disk objects (see also Brown et al. 2007). We varied these (and other) parameters 
in test simulations to determine the sensitivity of results to different assumptions. 

The SFD of the scattered disk was normalized as in Levison et al. (2008b). We as- 
sumed that there are presently 50,000 scattered disk objects with P > 50 km, and that this 
population decayed by a factor of 250 times since LHB (Tsiganis et al. 2005). We fixed 
/at = 1, because the dynamical depletion of CKBOs should be relatively minor during this 
stage. For R* = 37 km, this gives the initial mass ~0.2 Mg. Also, to promote collisional 
grinding, we used the weak disruption law with /g = 10. 

Figure [3] shows the SFD of CKBOs and scattered disk objects that were obtained for 
two different assumptions on the initial SFD of cold CKB. In Fig. |3K, we assumed that the 
initial SFD was steep down to P = 1 km. In Fig. [3b, we used the present SFD shape of cold 
CKB with R* = 37 km. 

In Fig. |3|l, the roll-over in the final SFD of CKBOs occurs at P ~ 5 km, which is 
significantly below the observed P* value. We experimented with a range of /at, /q, Vi, 
and Pj values, and various initial SFDs. These tests showed that the results illustrated in 
Fig. [3^ are representative. Specifically, if we start with the initial break at P < 10 km, 
the final break ends up being at P < 10 km as well. These results therefore suggest that 
the observed SFD roll-over of cold CKBOs at P* = 25-50 km cannot be produced by the 
collisional grinding after LHB. 
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On the other hand, if we start with the break at R* = 37 km (Fig. [Sb), the SFD of 
CKBOs remains nearly unchanged for R> 10 km. The total mass loss due to the coUisional 
grinding of small objects is only ~15% in the CKB and ~5% in the scattered disk. These 
results show that, while the SFD may have been shaped by collisional grinding in epochs 
prior to the LHB, it remained essentially constant since then. The binary fraction does not 
provide any interesting constraints on the collisional evolution of CKBOs after LHB because 
the vast majority of binaries with i? > 10 km survive. 

3.4. Pre-LHB Epoch 

Fraser (2009) and Benavidez & Campo Bagatin (2009) suggested that the observed 
SFD roll-over in the cold CKB was produced by collisional grinding during the ~600 My 
before the LHB. The collisional modeling of the pre-LHB phase is complicated by the fact 
that the state of the Kuiper belt before LHB is poorly understood. What we seem to infer 
from observations (see discussion in §1 and §3.2) is that the cold CKB probably formed in 
situ at 40-50 AU. The LHB modeling would then indicate that various dynamical processes 
probably removed ~90% of its mass during LHB (Morbidelli et al. 2008). Using these results 
as a guideline, we find that there is indeed a potential for the observed SFD roll-over being 
a fossil remnant of the collisional grinding of massive population before LHB. 

The impact speeds in the pre-LHB CKB depend on the dynamical state of the disk. 
In its present state, the impact speeds are relatively high (~1 km s~^) because orbits have 
relatively high eccentricities (e ~ 0-0.2) and inclinations {i ~ 0-30°). On the other hand, 
KBOs can only form if e and i were much lower (e.g., Kenyon et al. 2008 and the references 
therein). This implies that some dynamical processes must have excited e and i to their 
present values. Two main possibilities exist. Either (1) the cold CKB was dynamically 
excited by some primordial process that dates back to KBO formation, or (2) the excitation 
was produced by the LHB itself (see Morbidelli et al. 2008 for a discussion). Following 
Fraser (2009) we first study (1), in which case f i ~ 1 km s^^. Low Vi values implied by (2) 
will be discussed in §3.7. 

Fraser (2009) considered the initial SFDs that were steep (gi = 5) down to Ra, then 
became nearly flat {q2 = 0-2) down to Rf,, with the collisional equilibrium below R^ {q^, = 3.5). 
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These initial SFDs were motivated by the results of published simulations of the coagulation 
growth in KB, which tend to produce such distributions (Kenyon et al. 2008). Since it is not 
clear, however, whether KBOs formed by two-body coagulation in the first place (see Chiang 
& Youdin 2010 for a review), it is not guaranteed that these initial conditions actually apply. 
Nevertheless, we use Eraser's initial SFDs as a starting point and consider other options in 
§3.6. 

Figure m illustrates the Fraser's case with Ra = 2 km, Rt, = 0.5 km, gi = 5, g2 = 1 and 
ga = 3.5. As in Fraser (2009), we used Vi = 1 km s~^, fq = 3 (corresponding to Fraser's 
weak disruption law), and evolved the populations with Boulder over r = 600 Myr. Factor 
/at was varied to obtain different initial masses and thus different coUisional histories. Each 
of these cases would imply a different dynamical depletion factor during LHB. 

The best results were obtained with ~20-50, implying the initial mass between 7 
and 15 Me (Fig. Hi). With < 20 (<7 Me) and fjv > 50 (>15 Me), the SFD roll- 
over occurred at radii that were either too small and too large, respectively, compared to 
observations. These values are only soft limits, however, because the roll-over radius is also 
sensitive to the assumed Ra value. For example, smaller initial mass values would still be 
plausible if Ra = 3-10 km. In addition, slightly larger roll-over radii can be produced with 
fq = 10. We therefore find, in agreement with Fraser (2009), that a reasonably conservative 
lower limit on the initial disk's mass is ~1 Mg. 

The model SFDs discussed here share a common trait. While they bend to a shallow 
slope at R*, the slope below R* is never shallower than g ^ 3 and, even in the best cases such 
as the one shown in Fig. H^, just barely matches the observational constraint. Moreover, the 
shallow SFD segment below R* generally only extends down to R = 10-20 km, and steepens 
back to g ~ 5 for i? < 10 km (Fig. H^). This does not contradict the existing observations 
because very little is known about KBOs with i? < 10 km. We were unable obtain a case 
where the final SFD would be uniformly shallow below R* down to R < 10 km, except for 
very large and probably implausible initial masses. This is therefore clearly not the idealized 
case considered in PS05, where it was assumed that the slope below R* can be approximated 
by g ~ 3 to some very small (indefinite) R. 

Our simulations show that the disk mass is reduced by a factor of ~10 by coUisional 
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grinding. Thus, starting with ~1 Me, the remnant mass just before LHB would be ~0.1 Me- 
This is plausible because the LHB modeling in Morbidelli et al. (2008) showed that the 
population at 40-50 AU becomes dynamically depleted by a factor of ~10. The expected final 
mass of cold CKBOs from these order of magnitude considerations is therefore ~0.01 Me, 
which is in the right ballpark when compared to observations (e.g., Fuentes & Holman 2008, 
Brucker et al. 2009). 

The initial masses smaller than ~1 Me would imply Ra > 10 km, and indicate that 
the observed break at R* ~ 37 km was essentially in place before fragmentation processes 
had begun. The initial masses larger than ~10 Me are probably implausible, because these 
large masses would imply >1 Me mass before LHB and would require a dynamical depletion 
factor at LHB in excess of ~100. Such a large depletion factor is difficult to explain by LHB 
processes if the cold CKBOs formed in situ at 40-50 AU (Morbidelli et al. 2008). 

3.5. Binary Survival 

We now consider the binary survival. Figure shows the fraction of binaries surviving 
the pre-LHB epoch for /^r = 30 and Fraser's initial SFD shape. There is a clear trend with 
the physical size of binary components. Specifically, more than 50% of binaries with radii 
i? > 50 km survive, while the survival rate for R = 10-20 km is only ~0.5%. This trend 
is easy to understand because, according to Eqs. (!8l)- (|T0l) . smaller binary mass m implies 
larger orbital change. 

The trend is reversed for R < 10 km because of the lack of small impactors with 
R ~ 0.5 km that could unbind binaries with ~ 5 km (see Fig. H^). The behavior of 
the surviving binary fraction for i? < 10 km is sensitive to the initial SFD. For example, if 
Ra = Rb = '2' km, in which case the Dohnanyi's slope is directly attached to the steep SFD 
slope at large sizes, the survival rate of i? < 10 km binaries more monotonically drops with 
decreasing R (see §3.6). We will concentrate on binaries with i? > 10 km in the following 
discussion, because that's where things can be constrained by the existing observational 
data. 

Figure [5|l shows the physical parameters of known binaries in the cold CKB (Noll et al. 
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2008a, Grundy et al. 2009). The radii of primary components range from ~30 to ~100 km. 
The separations are between ~2 x 10^^ to ~0.2 Hill sphere. From Eq. ([1]), the survival rate 
should decrease with separation because 6v/v oc Vi^fa^ which is larger for larger a. Given 
the spread of separations in Fig. |5^, we therefore consider cases with a = 0.001, 0.01 and 
0.1 i?H- These initial values should cover the interesting range of separations. Considering 
these cases separately, rather than using some continuous initial distribution of a, is the right 
thing to do, because the real distribution of a produced by the formation process is not well 
understood. 

Figure ISJd shows the surviving binary fraction for these semimajor axis values. As 
expected, the wider binaries have lower survival rates than the tighter ones. For example, 
R = 30 km binaries with a = 0.1 -Rh would be reduced by a factor of ~100 in the Fraser's 
pre-LHB collisional model illustrated in Fig. H^, while R = 30 km binaries with a = 0.01 i?H 
would be reduced by a factor of ~10. These considerations provide an interesting test on 
the level of collisional grinding in the cold CKB, because the binary fraction could have been 
strongly reduced by collisions, especially for i? < 50 km. To pass this test, any plausible 
collisional scenario needs to match not only the SFD of cold CKBOs, including the roll-over 
at R* = 25-50 km, but also explain how the large fraction of cold CKBOs binaries survived, 
as indicated by observations (see §4). 



3.6. Sensitivity to Initial SFD 

Before we compare our results to observational data, we test the sensitivity of the binary 
survival to different assumptions. Figures |6] and [7] illustrate the dependence on the initial 
SFD. In Fig. [6l we choose to extend the steep slope (g = 5) from large sizes down to 
R = 2 km, and assume that there are no bodies initially with R < 2 km (case A). In Fig. 
[71 we impose a break from g = 5tog = 3at-R = 2 km (case B). In each case, /at, or 
equivalently the total initial mass, is set so that the collisional grinding leads to the SFD 
roll-over at i? = 25-50 km. The population is evolved over r = 600 My. 

The best results were obtained with Jn = 14 (4 Me initial mass) for case A and J'n = 30 
(26 Me) for case B. These initial masses grind down to 0.6 and 1 Me, respectively. While 
the final distribution obtained in case A matches constraints reasonably well (Fig. [6^, but 
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see discussion in §3.4), the distribution for R < 30 km in case B is always steeper than g = 3. 
Interestingly, it is difficult to produce a shallower slope in case B unless we use the initial 
masses in excess of 100 Me, which is clearly implausible, or fg » 10, which would conflict 
with the published results of impact simulations (see §3.1). 

It thus appears that a rather abrupt change in the initial SFD slope is needed to produce 
a shallow slope with q = 2-3 below R*. This can be easily understood. The objects with R 
near the initial slope change are long-lived, because they see a small number of impactors, if 
the transition is sharp. Thus, they can break larger bodies and create a sharp SFD roll-over 
at about ten times their radius. We find that it is easier to fit constraints if the initial break 
is placed at ~ 5 km, rather than at i? < 2 km, because i? ~ 5 km objects are long-lived 
in that case and can disrupt KBOs near R* with our disruption laws (/q = 1-10). 

Additional dependencies exist on the assumed SFDs of fragments produced by the cra- 
tering and catastrophic impacts. We find that catastrophic impacts tend to produce a sharp 
transition from steep to shallow slope near the largest object in the population that can be 
disrupted by them. The cratering impacts, on the other hand, tend to smooth this transition 
and produce more gentle waves in the SFD. This may explain some of the subtle differences 
between our simulations, which tend to produce gentle SFD waves, and those of Fraser 
(2009), which show stronger variations of slope for R R* (e.g., q < 2 below the break). 

The fraction of binaries surviving in case A is within a factor of ~2 to that we obtained 
with the Fraser's initial distribution. Interestingly, the binary survival rate is slightly larger 
in case B, which has larger mass than case A, and should thus lead to more perturbations on 
binary orbits. The larger binary survival rate in case B is related to the shorter lifetime of 
intermediate-size impactors (i? ~ 5 km) in case B, which are disrupted by smaller impactors 
from the Dohnanyi's tail. These small impactors are nearly absent in case A. 

3.7. Sensitivity to Impact Speeds 

So far we described tests with Vi ^ 1 km s~^. It is uncertain whether these assumed 
impact speeds should apply to the pre-LHB collisions, because the dynamical state of the 
Kuiper belt region at 40-50 AU could have been very different. Specifically, the dynamical 
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effects during LHB such as, e.g., passing resonances, could have excited the orbits of cold 
CKBOs. If so, the cold CKBOs could have had smaller eccentricities and inclinations than 
they have now, implying lower collision speeds before the LHB. Here we test cases with 
Vi < 1 km s~^. The results need to be considered with caution, because it is not clear 
whether our disruption laws (§2.1) are applicable with low impact speeds. 

Figure [H] illustrates the case with Vi = 300 m s"^. To create the roll-over at R* with this 
low collisional speed, we needed to assume fjy = 300, giving the initial mass of 88 Me for 
the case-A initial SFD. The population grinds down to 6.7 Me at r = 600 My, which would 
imply a dynamical depletion factor of >200. Collision speeds Vi < 300 m s~^ would require 
even larger initial masses and dynamical depletion factors that are clearly implausible. We 
therefore find that it is difficult to produce the SFD roll-over by collisional grinding with low 
collisional speeds (fj < 300 m s"^). 

The surviving binary fraction for Vi = 300 m s~^ (Fig. [8]d) shows the usual trend with 
R in that the binaries with physically smaller components survive at lower rate than the 
larger ones. Below R = 10-20 km, the surviving binary fraction drops to <10~'^. The results 
for i? > 20 km are similar to those obtained with f j ~ 1 km s~^ indicating that the survival 
of larger binaries is not overly sensitive to the assumption on f j. 

Additional tests show that it might be more plausible to create the roll-over at R* with 
low Vi, if ri ~ 5 km (initial 20 Me grinds to 5 Me) and/or for /q = 10 (20 Me grinds to 
2 Me). The results for binary survival are similar in these two cases. While ri = 5 km 
generates a slightly stronger gradient with R and minimum survivability for i? < 20 km, 
/q = 10 produces a shghtly softer gradient that is similar to that in Fig. [6)d. 

3.8. Changes of Binary Orbits 

So far we discussed the binary survival in different collisional scenarios. Here we describe 
the distribution of binary orbits of the surviving binaries. Figure [9^ shows the semimajor 
axis distribution of binary orbits produced in simulations with initial a = 0.01 Ru, and 
R = 20, 50 and 100 km. Figure [9b shows the same case for A^(a)da oc da/a for a < 0.1 -Rh, 
and N{a)da = for a > 0.1 Ru- These results were obtained for the Fraser's case (see 
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Fig.S]). 

While the binaries with R > 100 km tend to retain the shape of their original distribu- 
tion, the binary orbits for i? < 50 km become significantly modified. For example, Fig. |9^ 
shows that binaries with R = 50 km and a = 0.01 Ru become tighter or looser, producing 
a wide range of separations with a tail extending to a = 0.1 -Rh- The number of very wide 
binaries produced by collisions, however, is not expected to be large, because the tail of the 
distribution at a ~ 0.1 -Rh in Fig. ^ only represents a very small fraction. 

Figure IHb shows how the initial gradient of binary fraction with a can be modified by 
collisions. For example, the initial gradient oc for R = 20 km binaries, becomes oc a^^-^\ 
Thus, if the coUisional effects on these binary systems were important, it could be difficult 
to try to infer their primordial semimajor distribution from present observations. On the 
other hand, the semimajor axis distribution of binaries with R > 100 km should not have 
changed much since their formation. 

4. Comparison with Known Binaries 

All simulations that we conducted so far showed the following result. If the parameters 
were set up so that collisional grinding produced the SFD roll-over of cold CKB at R* = 25- 
50 km, the final binary fraction showed a strong gradient with radius. Therefore, such a 
gradient, if identified by observations of binary KBOs, would be a direct evidence for the 
extensive collisional grinding in the Kuiper belt. The absence of such a gradient, on the 
other hand, would indicate that the roll-over is more likely accretional. 

Figure [10] shows the fraction of binaries in the cold CKB that can be inferred from the 
existing observations (Noll et al. 2008a,b). The fraction shown here was roughly corrected 
for the main observational biases. For example, the secondaries can be more difficult to 
detect near small primaries due to their intrinsic faintness, and also because physically 
smaller components are expected to have smaller separations. We normalized each binary 
to Hy = 8.3 by scaling its component radii and separation by a factor. The binaries with 
normalized separations smaller than 0.032 arcsec, which was the smallest separation detected 
in the dataset, were removed. In total, only 5 out of 33 binaries were removed by this 
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procedure indicating that the observational bias is not overly important. 

The statistical errors shown in Fig. [10] are large, making it difficult to reach definitive 
conclusions. Still, some interesting features can be pointed out from data. First of all, the 
binary fraction for Hy < 5.5 {R > 100 km) seems to be lower than the one for Hy > 5.5. 
Since these large binaries would be relatively resistant to collisions during the solar system 
history, they paucity probably tells us something about the binary formation process. 

The fraction of binaries with Hy > 5.5 {R < 100 km) does not show any strong gradient 
with Hy (or R). Instead, the binary fraction is relatively constant and large down to the 
smallest surveyed objects {R ~ 30 km). In contrast, we found in §3 that the effects of 
coUisional grinding should deplete binaries with i? ~ 30 km, relative to those with R ~ 100 
km, by a factor of ~10. We therefore believe that the existing data are not suggestive of the 
kind of trends that we would expect to see in a population that experienced strong collisional 
grinding. This may indicate that the SFD roll-over in the cold CKB at R* = 25-50 km was 
not produced by disruptive collisions, but was instead already in place when KBOs were 
forming. 

Better observational statistics will be needed, especially for i? < 30 km, to test this 
preliminary conclusion. Additional caution needs to be exercised when comparing our colli- 
sional model results with observations, if the formation mechanism was capable of producing 
the binary fraction that strongly varied with R. For example, a relatively constant binary 
fraction could result from a combination of the formation and evolution effects, if the initial 
fraction was larger for smaller R and was modified by the collisional removal of small bina- 
ries. On the other hand, it is unlikely that the initial binary fraction was exactly 100%, as 
assumed here. If it were lower, our results could be used to place even harder limits on the 
extent of collisional grinding in KB. 

Another interesting feature in Fig. [10 is the dip at Hy = 6.8 [R = 60-70 km), where 
only 2 out of 31 surveyed cold CKBOs (~6%) turned out to be binary. Our collisional 
simulations were capable of producing such a dip (see, e.g. Fig. 4b), but at smaller radii 
{R = 10-20 km). We performed a Monte Carlo search in parameter space to identify cases 
that could produce the observed dip at R = 60-70 km. The best results were obtained 
with the initial break at i? > 10 km, Vi < 100 m s~^ and substantial initial mass. The low 
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coUisional speeds were required here so that the dip radius ended up to be only a factor of 
~2 larger than the final SFD roll-over radius. The SFD for it! > 10 km did not change much 
in these simulations implying that the SFD roll-over at R* — 25-50 km would have to be 
pretty much in place before the coUisional evolution started. 

5. Conclusions 

The work presented here shows how the Kuiper belt could have been affected by colli- 
sions. We found that extensive coUisional grinding, required if the SFD roll-over at R* = 25- 
50 km in cold CKB were coUisional, should imply a strong gradient of binary fraction as a 
function of R and separation. The current observational data do not show signs of such a 
gradient and instead suggest that small binaries (i? < 50 km) are at least as common as 
the large ones {R > 50 km). This may indicate that the SFD roll-over of cold CKBOs at 
R* — 25-50 km is not due to a prolonged phase of coUisional grinding in the Kuiper belt. 
Instead, the roll-over may be a fossil remnant of the KBO formation process. Future surveys 
of small binary KBOs will be able to test this conclusion. 
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Fig. 1. — Contrasting assumptions on the strength of KBOs. The plot shows the specific 
energy of catastrophic disruption, Q^, as a function of radius of the impacted body. The 
upper sohd curve shows as determined from hydrodynamic simulations of impacts on 
water ice (Benz & Asphaug 1999). The dashed curve shows Benz & Asphaug's ice 
divided by 10, which is the weakest disruption law considered in our work, and also roughly 
the smallest value in the gravity regime found in numerical simulations of impacts into 
porous rubble piles (Leinhardt & Stewart 2008, Stewart & Leinhardt 2008). The two solid 
straight lines show Qd that we computed from the strong (/3 = 3) and weak (/3 = 3/2) 
disruption laws of PS05. 
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Fig. 2. — Kuiper belt SFD changes for: (a) fq = 3, and (b) fq = 500. These two fragmen- 
tation laws roughly correspond to the strong and weak disruption laws considered by PS05 
= 3 and P — 3/2, respectively, in their notation). The upper thin black hues in both 
panels show the initial SFD (/^r = 1). The bottom bold black lines are the SFD produced by 
collisional grinding over r = 4.5 Gy. The bold grey lines denote constraints on the present 
SFD of the cold CKB: = 5 for R> 37 km according to Fraser et al. (2010) and g< = 2-3 
for R < 37 km for according to Fuentes et al. (2008, 2009). Since we plot the cumulative 
distributions A^(> R) here, the plotted slopes have indices equal to 4 for R > 37 km and 1-2 
for i? < 37 km. 
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Fig. 3. — SFDs of CKBOs (solid lines) and scattered disk objects (dashed lines) in our 
Boulder simulations of the post-LHB epoch for: (a) the initial SFD of CKBOs steep for 
i? > 1 km (g = 5) and shallow for i? < 1 km (g = 3), and (b) initial SFD of CKBOs with 
a break at R* — 37 km. In both cases, we assumed that the initially massive scattered disk 
dynamically decayed as in Tsiganis et al. (2005). The thin hues show the initial distributions. 
The bold lines are the final SFDs evolved over r — A Gy. We used fg — 10 and /n — 1. 
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Fig. 4. — (a) SFD evolution and (b) binary survival before LHB for the initial SFD shape 
used by Fraser (2009): Ra = 2 km, Rh = 0.5 km, gi = 5, g2 = 1 and = 3.5. We used 
fQ = 5 corresponding to the Fraser's weak disruption law. The thin black line in (a) shows 
the initial SFD. The middle bold black line in (a) shows the final state after r = 600 Myr 
of coUisional grinding. The bottom bold black line in (a) shows the SFD after a dynamical 
depletion factor /jv = 30 was applied to the final distribution (change indicated by an arrow) . 
The bold grey lines in (a) denote observational constraints on the present SFD in cold CKB. 
The bold grey line in (b) marks the binary fraction of 0.3. We assumed that the initial 
binary fraction was 1. The bold black line in (b) shows the expected binary fraction in this 
model. The binary separation was set to 0.01 Ru. 
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Fig. 5. — (a) Physical parameters of known binaries in the cold CKB (Noll et al. 2008a). 
Separations were computed from the apparent angular separation of the two components at 
discovery. Radius of the primary was estimated with j)y = 0.2. The horizontal dashed lines 
in (a) denote the separation values considered in our model, (b) The same as Fig. Hb, but 
for three different initial separations of binary components. As expected, the wide binaries 
have lower survival rates than the tight ones. 
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Fig. 6. — The same as Fig. Hlbut for the initial SFD that lacks objects with R < 2 km. We 
set Jn = 14 corresponding to the initial mass of ~4 Me. The initial population grinds down 
to ^0.6 Me in r = 600 Myr before LHB. 
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Fig. 7. — The same as Fig. H] but for the initial SFD with a break from q = 
a.t R = 2 km. We set /a? = 30 corresponding to the initial mass of ~26 Me. 
population grinds down to ^1 Me in r = 600 Myr before LHB. 
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Fig. 8. — The same as Fig. [6] but for Vi = 300 m s ^. We set = 300 corresponding to 
the initial mass of ~88 Me. The initial population grinds down to ~6.7 Me in r = 600 Myr 
before LHB. 
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Fig. 9. — The semimajor axis distribution of binary orbits for two different initial distribu- 
tions, (a) All binary orbits have a = 0.01 i?H initially (dashed line), (b) The number of 
binary orbits per interval da decreases as 1/a for a < 0.1 i?H (dashed line). We plot results 
for R = 20, 50 and 100 km for the Eraser's case illustrated in Fig. HI 
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Fig. 10. — The binary fraction in the cold CKB as a function of: (a) absolute primary 
magnitude Hy, and (b) primary radius. Data for 122 KBOs with i < 10°, including 33 
binaries, were taken from the observational surveys described in Noll et al. (2008a, b). The 
binary statistics for i < 5° is very similar to the one shown here. Crudely debiased data 
(see main text) were binned using AHy = 0.4. For Hy < 5.85, we combined all data into 
a single bin. The radii of objects in (b) were computed with pv — 0.2. The error bars and 
upper limits of non-detections show formal 1-sigma uncertainties. 



